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1 . 0 SUMMARY 


This report describes the hybrid-computer simulation of the General 
Electric OTW (Over The Wing) QCSEK (Quiet, Clean, Short-haul, Experimental 
Engine) control system. Tlie control system manipulates two variables: fuel 
flow and compressor stator angle. The system includes a full authority 
digital electronic control. The control system also contains a new Failure 
Identification and Corrective Action (FICA) scheme. The FICA logic is 
designed to recognize failures in the engine sensors of the digital control 
and to provide the corrective action for maintaining safe control of the 
engine. The digital control feeds a Fail-Fixed Servovalve (FFSV) contained 
in the hydromechanical fuel control portion of this system. The FFSV causes 
the fuel metering valve power piston to fail-in-position upon loss of elec- 
trical power or upon a hardover, electrical failure in either direction. 

The primary purpose of the simulation has been to develop this control 
system design with the following objectives: 

• Fast engine-thrust response for powered-lift operations. 

• Accurate steady-state and fast-response control of the engine where 
thrust variations are maintained within acceptable limits. 

Simulation results for throttle bursts from 62 to 100 percent net thrust 
predict that the experimental engine will meet the thrust requirement of 62 
to 95 percent thrust in one second. Operation of the digital control with 
the fail-fixed servovalve can be designed to meet the accuracy and fast- 
response objectives. Furthermore, results indicate that the digital control 
FICA strategy will continue to provide thrust management when several sensors 
have failed. 


2.0 INTRODUCTION 


The major purpose of the QCSEE program is to develop and demonstrate the 
technology required for propulsion systems of quiet, clean, economical, STOL 
(Short Takeoff /Landing) aircraft. One element of this program is to provide 
the digital electronic engine control technology necessary for an over-the- 
wing (OTW) turbofan engine application. 

Two key functions are to be performed by the digital electronic control. 
The first is to Improve the thrust response of a conventional turbofan to 
meet the stringent requirements of a STOL aircraft. The second function is 
to provide sensor failure Identification and corrective action (FICA) by 
using state-estimation techniques of modern control theory. In order to 
design this control and evaluate performance over the flight map, a hybrid 
simulation of the engine and control system was developed at the General 
Electric AEG (Aircraft Engine Group) Dynamic Analysis Simulation Center. 

This report contains a general description of the analytical models for 
the engine and control components which have been used in the simulation. 

The simulation was Implemented on the Simulation Center's two Electronic 
Associates, Inc. (EAI) Model 690 Hybrid Computing Systems. The engine and 
basic control system were constructed on one hybrid computing system to 
operate with or without FICA supervision from another hybrid computing sys- 
tem. This is possible because the FICA logic is designed to have no In- 
fluence on the control unless there is a sensor failure. 

Hybrid simulation data describing forward-thrust transients are pre- 
sented for sea level, static, standard-day conditions. Simulation results 
for a fall-fixed servovalve and a failure identification and corrective 
action control are also included. 
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3.0 OTW EXPERIMENTAL PROPULSION SYSTEM 


3.1 ENGINE 

The OTW experimental propulsion system, shown in Figure 1, features a 
high Mach (accelerating) inlet, a gear-driven fan, a composite fan-vane 
frame, an acoustically treated fan duct with an acoustic splitter ring, a 
variable-geometry confluent exhaust nozzle, an advanced core and low pressure 
turbine, an acoustically treated core exhaust, engine-mounted accessories, 
and a digital electronic control 5iystem. 

The fundamental design criterion which established the engine design 
approach was the fan engine cycle required to meet the noise objective. The 
fan and core exhaust pressure ratios were dictated by jet-flap noise con- 
straints and by the powered-lift requirements of an over-the-wing installa- 
tion. Analysis indicates that, when scaled in accordance with the specified 
ground rules, the engine will meet all of the program noise objectives. 

The fan is a low pressure ratio (1.35), low tip speed (350.5 m/sec - 
1150 ft/sec) configuration sized to provide 405.5 kg/sec (894 Ib/sec) of 
corrected airflow at takeoff power setting. It contains 28 titanium fan 
blades hav^ >g an aerodynamic contour capable of conversion to composite 
material in a flight system. The fan is driven by the low pressure turbine 
through a main reduction gear with a 2.0617 gear ratio. 

The "D" shaped, confluent exhaust nozzle incorporates side doors to vary 
the area from takeoff to cruise. The experimental engine will not Include a 
variable-exhaust-nozzle-actuation system; however, an actuation system for a 
flight-type engine can be extrapolated from current design technology. 


3.2 CONTROL SYSTEM 


The OTW engine control system manipulates fuel flovr and core compressor 
stator angle to achieve an excellent balance between thrust, fuel consump- 
tion, noise, exhaust pollution, and transient response. New technology 
features to be demonstrated on the OTW experimental engine control system 
are: a Full Authority Digital Electronic Control (FADEC) , a Failure Iden- 

tification and Corractive Action (FICA) scheme, and a Fail-Fixed Servovalve 
(FFSV). 

The OTW digital control is programmed to provide steady-state control, 
transient fuel schedules, core stator schedules, and limiting fuel-control 
functions. The transient fuel and core stator schedules on the UTW experi- 
mental engine were performed hydromechanlcTlly. 

Tha FICA strategy uses a form of Kalman filtering to estimate the con- 
trol system se.isor outputs and replaces sensor values with estimated values 
when deviations between the two exceed the expected maximum. 
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The fail-fixed servovalve Improves the reliability of the fuel delivery 
system. It causes the fuel metering valve power piston to fail-in-position 
upon loss of electrical power or upon a hard-over, electrical failure in 
either direction. 

The primary design requirements for the overall control system are; 

• Thrust control throughout the specified flight map with minimum 
pilot workload, 

• Fast thrust response - from 62 to 95% takeoff thrust in one second, 

• Specified noise and pollution goals. 

The general structure of the control system is shown in Figure 2. The 
digital electronic control is the heart of the system, manipulating fuel flow 
and core stator angle in response to a power command representative of an 
aircraft propulsion system computer. Digital control of fuel flow is main- 
tained through an electrohydraullc servovalve incorporated in a modified FlOl 
hydromechanical fuel control. In addition, the hydroroechaninal control 
provides the following: a backup core speed control, a positive metering- 

valve feedback to the digital control, a separate electrically operated 
emergency-shutdown function, and a pressurized fuel output from the fuel pump 
to power the conq>ressor stator servovalve and actuator assembly. 

Selection of a control mode which will accurately regulate thrust was 
based primarily on a tolerance analysis (Reference 1). The effects of typi- 
cal control and measurement tolerances, engine component tolerances, and 
engine component deterioration were determined for all potential control 
modes. The manipulated variables examined in this anlaysls were fuel flow 
and exhaust nozzle area, even though a variable exhaust nozzle is not in- 
cluded in the OTV experimental engine. Scheduling practicality, stability, 
and transient-response considerations were also factors in choosing the 
control mode. The control mode selected was one in which fuel flow controls 
corrected fan ^peed. 
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Figure 2. OT. Control Systfn, Schc»atl 











4.0 ANALYTIC.\L MODEL 


The analytical model represents the functional relations that exist 
between the variables of the OTW QCSEE and control system. The engine por- 
tion of the simulation Is based on the steady-state conservation equations, 
of the cycle deck, with the addition of significant dynamics. The model has 
detailed representations of the primary control system components: the fuel 

control and the compressor stator control. 


4.1 ENGINE MODEL 


4.1.1 General 

The engine model Is derived from the cycle deck which was used to gen- 
erate the QCSEE preliminary technical requirements. Significant Items repre- 
sented are rotor dynamics, heat soak, and combustion delay. The effects of 
heat soak, combustion delay, and compressor stator error are based upon past 
experience with the FlOl core engine. No distortion or relngestlon effects 
are Included; volume dynamic effects In the frequency range of consideration 
are Insignificant and have been omitted. 

The philosophy behind the model construction was to develop mathematical 
representations, for the components, to duplicate the thermodynamics used in the 
cycle deck. Memory limitations in the hybrid computer require linear table- 
interpolations; whereas, some cycle deck maps use polynomial interpolation. 

Due to the computation time and memory involved, the Iteration procedure 
of the cycle deck was considered undesirable. To keep cycle time at a mini- 
mum on the hybrid (maximum sampling rate), no actual iterative process in- 
volving tolerances is used; instead, all variables requiring an Iteration are 
computed using values based on the previous time sample. This results in two 
typer of iterative process: one internal to the component and one between 

jmponents. Since the latter cannot be evaluated without the components, the 
first step was to develop the component representations and internal "itera- 
tions," The choice of input and output variables was based on past experi- 
ence. With the exception of the fan and compressor map input variables, all 
of the component internal iterations can be considered as direct-substitution 
iterations under steady-state conditions. The component representations were 
considered adequate if, with inputs held constant, the iteration variables 
were stable and the outputs within two percent of the cycle deck. 

Based on past experience, it was recognized that some Iteration vari- 
ables might require filtering due to the iterative process between compo- 
nents. For this reason, a digital filter routine was constructed to produce 
a lag in the input variable based on an Euler integration using the average 
sampling rate. This changes the effective iteration process from direct 
substitution to weighted averaging under steady-state conditions. Variables 
requiring this treatment are fan front-face total flow used in the Installed 
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inlet, combustor discharge pressure In the HP (high pressure) turbine, and LP 
(low pressure) turbine vane inlet total pressure. Based on 16.67 samples/ 
sec and a 7.0:1 time scale, the respective time constant are 30, 12, and 48 
milliseconds. 

Figure 3 illustrates the information flow between components. 


4.1.2 Inlet 


The inlet simulation receives flight altitude, (ALT), Mach number (XM), 
and temperature increment (DTAMB) from the standard day ambient temperature 
as determined by the desired flight condition. Fan front-face total flow 
(W2A) is received from the fan when the inlet is in the Installed mode. Out- 
puts are free-stream pressure (PO) and temperature (TO), inlet throat static 
tPSll) and totaJ pressure (Pll), and fan front-face total temperature (Til) 
and pressure (P12). A constant specific heat ratio of 1.4 is used in the 
inlet simulation. No distortion or relngestlon effects are considered. The 
only dynamic effect Is a lag of 30 milliseconds in fan front-face total flow. 
This lag is necessary to provide model Iteration stability between components 
in the installed mode which produces a total pressure drop between the inlet 
throat and fan front face. 


4.1.3 Fan Tip 

The fan tip simulation receives front-fare total temperature (Til) and 
pressure (PI 2) f’om the inlet, physical speed from the LP rotor (XNL) , and 
tip discharge total pressure (P13) from the bypass duct. Outputs are tip 
discharge specific heat (CP13), tip inlet (H12) and discharge enthalpies, 
percent fan corrected speed (I’CNLR), tip inlet entropy function (PHI12), tip 
inlet gas constant (R12), tip disclmrge total temperature (T13), and front 
face total flow (W2A). No dynamics are used. To conserve memory, the poly- 
nomial, fan-tip maps used in the cycle deck have been converted to use linear 
interpolation. The general-parameter map input is an internal iteration 
variable; the value depends on the previous-pass value, discharge total 
pressure input, and the discharge total pressure computed from the map on the 
previous pass. 


4.1.4 Fan Hub 

The fan hub simulation receives front face total temperature (T12) and 
pressure (PI 2) from the inlet; fan tip inlet (H12) and discharge enthalpies 
(H13), inlet entropy function (PHI12), inlet gas constant (R12), percent fan 
corrected speed (PCNUl), and front-face total flow (W2A) from die fan tip; 
and compressor inlet airflow (W25) from the compressor. Outputs are fan hub 
(PW2) and tip (PW12) power, fan tip inlet airflow (W12), and compressor inlet 
total temperature (T25) and pressure (P25). No dynamics are Involved with 
the fan huh simulation. Reynolds number effects are the same as in the cycle 
deck. Ue.causc of memory limitations, the polynomial, fan-hub maps used in 
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the cycle deck have been converted to use linear interpolation in the hybrid 
tBodel. The general-parameter map input is an Internal iteration variable 
whose value depends on the past value and the difference between [hub inlet 
airflow (W12)]/[hub pressure ratio (P23Q2)] input and the previous map value; 
this uses the cycle deck assumption that all hub flow goes into the compres- 
sor. The total pressure drop between the fan hub discharge and the compres- 
sor inlet is calculated as in the cycle deck. 


4.1.5 Bypass Duct 

The bypass duct simulation receives fan tip discharge specific heat 
(CF13) and total temperature (T13)> and inlet gas constant (R12) from the fan 
tip; fan tip discharge flow (W12) from the fan hub; and bypass duct discharge 
total pressure (P16) from the mxler. Output is bypass duct inlet total 
pressure (P13). No dynamics are Involved in the simulation. Since the cycle 
deck variation in the fan tip discharge and the bypass uct inlet total 
pressure ratio is less than 0.00003, this pressure ratio Is assumed constant. 
The total pressure drop between the bypass duct inlet and discharge is cal- 
culated as in the cycle deck. 


4.1.6 HP Compressor 

The HP compressor simulation receives compressor inlet total pressure 
(P25) and temperature (T25) from the fan hub, compressor discharge total 
pressure (P3) from the combustor, HP compressor physical speed (XNH) from the 
HP rotor, and off-stator error (DLBI CA) from the control. Produced are 
compressor discharge (W3) and inlet (W25) airflows, combustor inlet gas flow 
(W36), all cooling flows and enthalpies, compressor discharge (H3) and com- 
bustor inlet enthalpies (H36), compressor discharge specific heat (CP3), 
compressor discharge (T3) and combustor inlet total temperature (T36), and 
compressor power (PW25) . The simulation uses the same base compressor maps 
used on the UTW simulation. Changes were made to agree with cycle deck 
modifiers and to be compatible with the new hybrid software. Reynolds number 
effects are assumed negligible and are omitted. The UTW simulation of off- 
schedule stator effects on airflow and efficiency was modified to provide an 
accurate representation for the large, off-schedule positions associated with 
reset of stators for OTW approach power transients. This modification was 
based on recent compressor test data. 

The only dynamics are due to heat soak producing an effective lag in 
compressor discharge total temperature. 


4.1.7 Combustor 


The combustor simulation receives compressor discharge specific heat 
(CP3), enthalpy (H3), total temperature (T3) , and flow (W3) from the compres- 
sor; combustor inlet enthalpy (H36), total temperature, (T36) and flow (W36) 
from the compressor; compressor inlet airflow (W25) from the compressor; fuel 
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flow (WFM) from the control; and combustor discharge total pressure (P4) from 
the HP turbine. Outputs produced are HP turbine rotor Inlet fuel/alr ratio 
(FARs'), enthalpy (HAl), total temperature (T41), airflow (WA41), gas flow 
(W41) and entropy function (PH141); and compressor discharge static (PS3) and 
total pressures (1’3). The total pressure drop between the compressor dis- 
charge (P3) and the combustor dlsharge (P4) is represented as in the cycle 
deck. The efficiency map Is Identical to the cycle deck. Two dynamic effects 
are included in the combustor; one is a 25-millisecond combustor delay be- 
tween the fuel flow produced by the control and the fuel flow that Is burned 
in the engine. The second Is a heat soak calculation which produces an 
effective lag In HP turbine rotor inlet total temperature (T41). Both 
effects are based on previous experience with the FlOl core engine. Mixing 
of primary gas flow and a cooling flow is accomplished for the HP turbine 
rotor inlet using steady-state conservation of energy. There Is no effect on 
total pressure. 


4.1.8 HP Turbine 


The HP turbine receives HP turbine rotor inlet fuel/alr ratio (FAR41) , 
entropy function (PH141), total temperature (T41), enthalpy (H41), gas flow 
(W41), and airflow (WA41) from the combustor; compressor discharge total 
temperature (T3) and percent compressor physical speed (XNH) from the com- 
pressor; cooling flows and enthalpies for the HP turbine discharge (WCL42), 

LP turbine vane inlet (WCL48), and I.P turbine rotor inlet (WCL49) from the 
compressor; LP turbine vane inlet total pressure (P48) from the LP turbine; 
and HP compressor physical speed (XNH) from the HP rotor. Outputs are com- 
bustor discharge i Jtal pressure (P4); HP turbine power (PW41); U turbine 
rotor inlet fuel/alr ratio (FAR49) , enthalpy (H49), entropy function (PHI49), 
total temperature (T49), gas flow (W49), and airflow (WA49); Id* turbine vane 
inlet total temperature (T48) and gas flow (W48). The component maps used are 
identical to those of the cycle deck. Mixing of primary gas flow and cooling 
flow is accomplished using steady-state conservation of energy at the HP 
turbine discharge, the LP turbine vane inlet, and the U turbine rotor inlet. 
There is no effect on total pressure. The only dynamic effect is a 12- 
milllsecond lag on combustor discharge total pressure, which is necessary for 
model iteration stability between components. 


4.1.9 L P Turbine 

The LP turbine receives LP turbine vane inlet total temperature (T48) 
and gas flow (W48) from the HP turbine; I.P turbine rotor inlet fuel/air ratio 
(FAR49), enthalpy (H49), entropy function (PH149), total temperature (T49), 
gas flow (W49), and airflow (WA49) from the HP turbine; LP turbine discharge 
cooling flow (WCL5) and enthalpy (HCL5) from the compressor; percent fan 
corrected speed (PCNLR) from the fan tip; fan physical speed (XNL) from the 
LP rotor; and LP turbine discharge total pressure (P5) from the mixer. Pro- 
duced are LP turbine power (PW49); LP turbine discharge fuel /air ratio (FARS), 
enthalpy (H5), gas flow (W5), and airflow (WAS); and U turbine vane inlet 
total pressure (P48). Ihe component maps used are identical to those of the 


cycle deck. Mixing of primary gas flow and cooling flow at the LP turbine 
discharge la accomplished using steady-state conservation of energy. There 
Is no effect on total pressure. As In the cycle deck* If turbine rotor and 
vane inlet total pressures are equal. The only dynamic effect is a 48- 
mllllsecond lag on LP turbine rotor inlet total pressure which Is necessary 
for model Iteration stability between components. 


4.1.10 Mixer 


The mixer' simulation receives fan tip discharge specific heat (CP13)* 
enthalpy (H13), and total temperature (T13) from the fan tip; fan tip inlet 
gas constant (R12) from the fan tip; fan tip Inlet airflow (W12) from the fan 
hub; LP turbine discharge fuel/air ratio (FAR5), enthalpy (H5), gas flow 
(W5)» and airflow (WAS) from the LP turbine; LP turbine frame discharge 
cooling flow (WCL5S) and enthalpy (HCLS5) from the compressor; and mixer 
discharge total pressure from the exhaust nozzle. Outputs are mixer dis- 
charge fuel/alr ratio (FAR6), specific heat ratio (GAM6), enthalpy (H6), 
total temperature (T55)* and gas flow (W6); bypass duct discharge total (P16) 
and static pressures (PS16) ; LP turbine discharge total pressure (PS) ; and 
LP turbine frame discharge total temperature (TSS). Mixing of primary gas 
flow and cooling flow at the LP turbine frame discharge is accomplished using 
steady-state conservation of energy. LP turbine swirl effect is represented 
as a constant (LP turbine frame discharge total pressure) /(LP turbine dis- 
charge total pressure). Mixing of the bypass and core flow is performed 
using steady-state conservation of mass, energy, and momentum. This also 
uses the cycle deck assumption that core duct discharge static pressure 
equals bypass duct discharge static pressure. No dynamics are involved. 

Total pressure drops from mixer discharge to mixing plane and core duct 
discharge to LP turbine frame discharge are represented as in the cycle deck. 


4.1.11 Exhaust Nozzle 

The exhaust nozzle simulation receives mixer discharge specific heat 
ratio (GAM6), total temperature (T6), and t;as flow (W6) from the mixer; 
exhaust nozzle throat actual area from the control (A8); and free-stream 
pressure (PO) from the inlet. Produced are exhaust nozzle throat effective 
area (AES), static pressure (PSS), and velocity head; and mixer discharge 
total pressure (P6). No dynamics are used. Nozzle flow coefficient maps are 
identical with those of the cycle ueck. 


4.1.12 Rotor Dynamics 

Rotor speeds are computed using conservation of angular momentum. The 
LP rotor receives fan tip (PW12) and hub powers (PW2) from the fan hub, LP 
turbine power from the LP turbine (PW49) , and a power-loss term that is a func- 
tion, derived empirically, based on cycle deck data. All moments of inertia 
for the LP rotor have been reflected to the fan side of the gearbox. Using 
this inertia, a dynamic value for fan physical speed (XNL) Is calculated. 
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For the HP rotor, a dynamic HP compressor pliysical speed (XNH) Is cal- 
culated using the core inertia and powers from the compressor (PN25), HP 
turbine (PW41), HP rotor loss calculation (PWPXH), and desired customer power 
takeoff. The HP rotor power-loss calculation is based on an empirical fit of 
cycle deck data. 


A. 2 CONTROL SYSTEM MODEL 

4.2.1 General 

In order to develop the control independent of FICA logic, the control 
system is simulated on the same EAT hybrid computing system as the engine. 

For the purpose of discussion, the control system model is divided into two 
sections: that which nutnlpulates core engine fuel flow and that which posi- 

tions the HP compressor stators. Both are primarily contained on the analog 
portion of the hybrid computer. As a consequence, the simulation does not 
account for round-off errors associated with the 12-bit words in the digital 
control computer. It has been assumed that the software program for the 
digital control computer will be designed to prevent deteriorating effects of 
round-off errors on control performance. 

Both simulations Include an analog approximation for the effect of the 
digital computer time delay. Hsiug the estimated cycle time of 10 milli- 
seconds for the digital control, the effect ol the digltal-to-analog conver- 
sion may be approximated by a 5-mll 1 isecond delay in the frequency domain 
(Reference 2). Therefore, the analog ot the digital control computer time 
delay is represented by a total of 15 mi 1 1 iseconds . Duo to equipment limita- 
tions, this time delay is further approximated by a first-order lag with a 
15-mlllisecond constant. This lag provides a reasonable approximation tor 
phase shift, in the range of the control loop crossover frequencies, which is 
important for stability considerations. The lag does provide some effect of 
time delay during large transients, but is not precise. In the simulations, 
the lag is located just prior to the torque motor driver amplifier simula- 
tions. 


4.2.2 Fuel (\^ntrol 

A schematic of the engine tuel control moilel is shown in ligure 4. Ihe 
model Includes representations for the digital electronic and the hydro- 
mechanical portions of the fuel iH>ntrol. The model for the digital electro- 
nic portion is based on the detailed block diagrams contained in Appendix A. 
Appendix A also contains the sped 1 icat ions lor the gains, time constants, 
and limits which were developed from the hybrid simulations of the engine and 
control system. Currently, these diagrams and specif icat ions are being used 
to program the digltiil luel I'ontrol lor t hi' 1 irst build ol the OlV experi- 
mental engine. It should be noted that the actual hybrid simulatiiin of the 
digital fuel control is not a one-to-one representation ol all details in the 
Appendix A block diagrams and spec ll icat ions . For example, it does not in- 
clude the polynomial equation curve-fits tor tlie accelerat ioji tuel schedule; 
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Figure 4. OTW Core Engine Fuel Control Model Schematic. 















the hybrid simulation contains the base functions for this schedule, from 
which the polynomial equation fits were derived. In summary, the hybrid 
simulation includes those details which were necessary to develop and evalu- 
ate steady-state and transient performance of the engine and digital fuel 
control . 

The model for the hydromechanical portion of the fuel control is defined 
by Figure 5. The hydromechanlcaJ portion is a modified FlOl fuel control. 

The model Includes representations for the backup core-speed-limit control, 
two servovalve designs, the metering valve newer piston, the metering valve 
port shape, and an effective 0.02-second leg for the delta pressure regulator 
that functions to maintain a constant pressure drop across the metering valve 
port. The two servovalve designs contained in the model are the fail-fixed 
servovalve and the linear servovalve normally used in the FlOl fuel control. 
The details of the fail-fixed servovalve are discussed in the next section of 
this report. The following traasTer function for the linear servovalve is 
used in the model; 

3 

A Servovalve Flow _ 0. 0191 in. /sec 

A Torque Motor Current 0.03 S + 1 mA 


where torque motor current is limited to the -80 to +80 mA range. 

As Indicated by Figure 4, several control modes are accommodated by the 
fuel control. In the power setting range from approach to takeoff, the pri- 
mary control mode is the manipulation of fuel flow (WF) to control corrected 
fan speed, which is a variable related to thrust. Tlie control mode is 
changed to either maximum core speed, maximum core turbine inlet temperature 
(T41C), accel schedule (WF/PS3\ or the decel schedule (WF/PS3) when the 
engine tries to operate beyond the scheduled limit of one of these variables. 
In addition, an inlet duct Mach number limit controller would be used in the 
digital fuel control design for a flight-type OTW engine. It is not Included 
for the experimental engine in order to reduce digital memory requirements. 
The actual hybrid simulations of the fuel control contains this controller 
for analytical purposes; however, it is not defined in the Appendix A block 
diagram and specifications. 

When the engine operates within tlie above scheduled limits, a fan speed 
error term is generated based on sensed fan speed and an uncorrected fan 
speed from the corrected fan speed schedule. The fan speed error and the 
lagged rate of change of sensed, metering-valve position are used to deter- 
mine the torque motor current output (IWT) from the digital control. This 
current positions an electrohydromechanical servovalve which Is ported to the 
metering-valve, power piston. The magnitude and polarity of the current 
determine the rate and direction of metering-valve, power-piston position 
(XMV). The metering-valve area is proportional to the square of the power- 
piston position. A pressure regulator maintains a constant pressure drop 
across the metering valve; tlms, the metered fuel flow to the engine is 
proportional to the square of the power-piston position. 
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For small perturbations, the transfer function from fan speed error to 
fuel flow Is approximated by: 


A Fuel Flow ^ K (0.3 S+1) 

A Fan Speed rpm S 


The gain "K" varies with respect to the square-law shape of the mter- 
ing-valve area. The lead-time constant in this transfer function is due to, 
and thus equal to, the time constant of the lag on the rate of change of 
sensed, metering-valve position. The valve of the lead-time constant is 
chosen to compensate for the engine lag from fuel flow to fan speed. In 
summary, the fan-speed/fuel-flow control is an integrating (i.e.. Type 1) 
control with lead compensation and is designed to provide accurate control of 
corrected fan speed. 

The dynamic designs of the limit controls for core speed, core turbine 
inlet temperature, and inlet Mach number are similar to those of the fan speed 
control; each is an integrating control with lead compensation. To com- 
pensate for engine and sensor lags, lead-time constants are 0.5 seconds for 
the core speed and inlet Mach number controls, and 0.1 second for the core 
turbine inlet temperature control. 

During engine accelerations, fuel flow is limited by the WF/PS3 accel 
schedule which is a function of HP compressor (core) speed and HP conq>ressor 
inlet temperature (T25). HP compressor inlet temperature is calculated based 
on sensed fan speed and sensed fan inlet temperature (T12). The accel sched- 
ule is designed to provide adequate HP compressor stall margin and to be a 
backup for the core turbine inlet temperature control. During engine decel- 
erations, fuel flow is limited by a constant WF/PS3 decel schedule to prevent 
combustor blow out. 


4.2.3 Fall-Fixed Servovalve 


The model of the fall-fixed servovalve used in the simulation studies 
included representations for two types of pulse-width-raodulated, torque- 
raotor-drlver amplifiers. The unipolar-pulse-drlver amplifier uses current 
pulses of only one polarity chosen by the sign of the digital word. The 
bipolar-pulse-driver amplifier uses a train of positive and negative pulses; 
pulse-width is determined by the sign and magnitude of the digital word. 

Both operate at a frequency of 500 hertz. Operation of the servovalve with 
the unipolar amplifier was felt to have better resolution from percent 
digital-word-to-flow, but poorer null shift and dead-zone characteristics 
than operation with the bipolar amplifier. Figure 6 Illustrates the differ- 
ence between servovalve operation with the unipolar and the bipolar ampli- 
fiers. 

Since the primary concern was the effect of the dead-zone and null 
shift, the dlgltal-word~to-;<ervovalve flow characteristics were linearized 
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Figure 6. P«ll-Flxed Servovalve Characterlatlca; Digital Anpllfler 
and Servovalve, 500 Hs, Pulae-Wldth Modulation. 
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and United aa ahovn in Figure 7. Tolerancea on the dead-tone size and 
Mount of null ahift were estinated baaed on paat experience with ainilar 
valvea. The pulae-driver ai^lifiera and the aervovalve were aaautMMl to 
reapond like firat-order laga with tine conatanta of 0.01 and 0.03 acconda 
reapectively. The arae tine conatanta were aaauaed for the torque-notor- 
driver anplifier and linear aervovalve diacuaaed in the preceding aection of 
thia report. 


4.2.4 Core Conpreaaor Stator Control 

The model for the core conpreaaor atator control la deacribed in Figure 
8. It Includea repreaentationa for the digital electronic and hydronechanl- 
cal portiona of the atator control. The model for the digital electronic 
portion la baaed on the detailed block diagram in Appendix B (Mote: See 

Appendix A for procedure to calculate T25 and corrected core apeed which ara 
uaed in both t\M fuel and the atator controla). Appendix B alao containa the 
apeclflcationa for the gaina. tine conatanta. achedulea. and llmlta which 
were developed from dealgn atudiea on the engine and control ayatera. Cur- 
rently, theae diagrama and apeclflcationa are being uaed to program the 
digital core conpreaaor atator control for the first build of the OTU experi- 
mental engine. It ahould be noted that the actual hybrid alraulation of the 
digital atator control la not a one-to-one repreaentatlon of all detaila in 
the block diagram and apeclflcationa. For exaaple, it doea not include the 
logic asaociated with detecting a failure in one of the poaition-feedback 
aenaora. However, the hybrid almilation doea include thoae detaila which are 
neccaaary for proper evaluation of ateady-atate and tranalent performance. 

The model for the hydromechanics 1 portion of the core atator control is 
defined in Figure 9. The velocity of the stator actuators, in units of cm 
(inches) ZBETA feedback stroke to the LVOT (linearly variable, differential 
transformer), is represented by the product of servovalve current (IB) times 
a linear function of fuel flow. This linear function is used to approximate 
the effect of the varying fuel-pinap-discharge pressure which is the supply 
source for the core stator servovalve. For example, an increase in fuel flow 
causes an increase in pinop-diacharge pressure and, thus, an increase in gain 
from servovalve current to flow output. The net result is a larger velocity 
gain from servovalve current to actuator position. The model also includes 
the nonlinear linkage relation between actuator position and actual core 
stator angle. 

T!te combined operation of the digital and hydromechanical portions of 
the core stator control in Figure 8 is as follows. When the digital control 
logic is set for the "reset off" mode, the output of the reset control is 
zero. Therefore, the position of the stator vanes is scheduled as a function 
of core corrected speed. The error between the scheduled and sensed stator 
positions determines the magnitude of torque motor current (IB) from the 
digital control. Tltis current determines the servovalve flow to the two ram 
actuators which position the stator ansle. 


When the digital control logic is set for the "reset on" mode, the reset 
control adds a scheduled amount of degrees to the nominal core stator sched- 
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ule. The net effect is to set the stators at a more closed position. For 
the OTW experimental engine, the core stator reset schedule is an adjustable 
function of power setting. The rate of adding reset is limited by the reset 
control. The rate of removing reset (i.e. , opening stators) is also limited. 
The details of the reset control are described by the block diagram in 
Appendix B. 


4.3 FAILURE IDENTIFICATION AND CORRECTIVE ACTION SIMULATION 

Failure Identification and Corrective Action (FICA) is an integral part 
of the digital control. The FICA will replace the output of a sensor in the 
event of failure. The FICA has no effect on unfailed-sensor output, so it 
can be separate from the rost of the control system. The inputs to the FICA 
are the sensors and the control outputs. The outputs from the FICA are the 
estimated values of the failed sensors and the values of the unfailed sen- 
sors. 


The FICA strategy has been developed and has first application in the 
OTW QCSEE full authority digital control. It is based on an extended Kalman 
filter incorporating a nonlinear model of the engine to provide a best esti- 
mate of the state and expected sensor outputs of the engine controls. 

Since the FICA acts only in the event of a sensor failure, it has no 
effect on the normal control action. Ute logic, schedules, and dynamics of 
the control were designed independently of the FICA, and the FICA was de- 
signed with the control in place. The simulation of the FICA was run on the 
Number 2 hybrid cumputtr while the engine and control of the QCSEE was simu- 
lated on the Number I hybrid computer. The connections between the two 
computers were analog trunk lines. In adding the FICA, the simulated sensor 
outputs were reconnected to the FICA, and the FICA supplied the sensor out- 
puts to the controls simulation. For unfailed-sensor operation the FICA does 
not alter the sensor output signals. 

The control inputs and sensed output variables for the QCSEE are shown 
In Figure 10. The inputs are current to the fuel flow valve, IWF, and cur- 
rent to the compressor stator blade torque motor, IB. The actual fuel flow 
and actual compressor stator angle are sensed output variables. Also sensed 
are the high and low rotor speeds, compressor discharge pressure and tempera- 
ture, and turbine discharge temperatures for a total of seven sensed output 
variables. In addition, the three environmental variables P2, T2, and PO are 
sensed and used for model Inputs. 

An overall block diagram of the failure detection and correction strat- 
egy is shown in Figure 11. The strategy consists of three parts: a non- 

linear model of the engine, decision logic to determine when a failure has 
occurred and to take corrective action when u failure is detected, and a 
feedback-gain matrix to update the model and keep it in close agreement with 
the actual engine. The engine control logic la unaffected by the failure 
detection and correction logic. The same control signals are applied to the 
model as well as to the engine. The outputs of the models, which are the 
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expected outputs of the engine, are compared with the actual outputs of the 
engine. If the difference between the actual and expected outputs Is too 
large. It Is assumed that the sensor has failed. In that case, the sensor 
output Is replaced by the expected output from the model for control pur- 
poses. The difference between the actual and expected output for an unf ailed 
sensor is fed-back through a gain matrix to correct the model. 

Obviously, this strategy relies heavily on the nonlinear model. This 
will be discussed In more detail below. The theory behind the feedback gain 
and decision logic Is also given below. Since the strategy Is to be imple- 
mented with an on-englne digital computer, the feedback matrix could not be 
computed in an on-line basis. This matrix is constant and Is not changed for 
any combination of sensor failures or engine-operating conditions, the 
strategy has been extensively tested using an accurate, hybrid simulation of 
the engine. Some results of these tests are shown under Simulation Results, 
Section 6. 


4.3.1 Model 


The model is designed to have good dynamic and static accuracy and yet 
be simple enough to be solved by an on-englne, digital computer. The high 
accuracy is necessary to allow good, tight control even though a sensor has 
failed. The engine model for the subsonic transport Is designed to be 
accurate for a power range from flight Idle to maximum, over an inlet tem- 
perature range from about -54 to 54® C (-65 to 130® F) and an Inlet pressure 
range from about 20 to 132 kN/m^ (0.2 to 1.3 atm) to accommodate a sensor 
failure at any condition In flight. 

The engine model has the same Inputs and outputs as the actual engine, 
as indicated by Figure 10, plus the inlet air pressure and temperature and 
the ambient air pressure that determine the external operating environment. 

The form of the model follows the engine cycle schematically shown in 
Figure 10. It consists of a low pressure ratio fan driven by the low pres- 
sure turbine that gets the exhaust gasses from the high pressure turbine. A 
small fraction of the fan discharge air enters the high pressure compressor 
that is driven by the high pressure turbine. Fuel is burned in a combustor, 
using air from the high pressure compressor, exhausting into the high pres- 
sure turbine. Most of the fan discharge air bypasses the compressor, com- 
bustor, and turbines and is mixed with the low pressure turbine exhaust 
gasses. The mixed gasses are accelerated in the Jet nozzle and exhausted. 

The model must account for the steady-state and transient performance of 
the engine. This is done in a simplified accounting for the mass flows, 
pressures, and temperatures in the engine and for the power in the fan, 
compressor, and turbines. The representation of the fans, compressors, and 
turbines is designed to approximate the operating line for steady-state and 
the off-operating line for transient response. These representations are in 
the form of polynomials and tables, and the choice in each case was made to 
get the simplest calculation with the desired accuracy. Rather than present 
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the detailed polynomials and tables, which will vary from engine to engine, 
the basic engine relationships will be presented In functional form Indicat- 
ing the Interrelationships between variables. The block diagram shown In 
Figure 12 shows how these relationships follow from the engine cycle of 
Figure 10. The nomenclature Is: W - gas flow, P - stagnation pressure, T - 
stagnation temperature, and the following characters, numbers, or subscripts 
provide Identification (see Figure 12 and Appendix C) . Since there Is a 
significant lag In the temperature sensors, the output of these sensors Is a 
state and Is Indicated by a terminal S In the nomenclature. 

For the fan, the airflow Is of the form 


W2 - (P2/*^) fj^ (NL/^^2, P15/P2) (1) 

and the temperature rise across the fan is 

T15 - T2 X fj,^ (NL/»^, P15/P2) (2) 

where NL Is the fan speed, rpm. 

The core inlet pressure and temperature are affected by the flow split 
and core Inlet ducts. 

T21 - fj^ (T15, T2) (3) 

P21 - P15 X fjp (NL/r^) (4) 

The high pressure compressor airflow is 

W21 - (P21/»'^n) X (NH/*^, PS3/P21, Z3) (5) 

where NH is the high pressure rotor speed, Z(^ is the stator position, and the 
temperature rise is: 

T3 = T21 X (NH/*'t2. PS3/P21 (6) 

The combustor Inlet airflow Is 

W3 « C„ X W21 (7) 

w 

where C,, Is a constant. 

W 

The combustor teroperatuie rise depends on fuel/alr ratio 


Actuator 



Figure 12. Block Diagram of Engine Model Computation. 















T4 


T3 + (WF/W3) 


( 8 ) 


where WP is the fuel flow rate and discharge airflow is: 


W4 - W3 + WF (9) 

The high pressure turbine inlet pressure is based on a fixed-area nozzle. 

P4 - C^„ X W4 X (10) 

TW 

The combustor inlet static pressure is related to the turbine inlet 
pressure. 


PS3 - C_, X P4 (H) 

oL 

The fuel flow is a function of the metering valve position. 

WF - (Z^) (12) 

The bypass duct airflow is the fan flow less the compressor airflow. 
W15 = W2 - W21 (13) 


The turbine discharge pressure is computed from the bypass duct pressure 
loss. At the mixing plane between the core gas flow and bypass airflow, the 
approximation that the static pressures of the two streams are equal is used. 


PS56 

- P15 - f„, (W15, P15, T15) 

dL 

(14) 

W56 

- W21 + WF 

(15) 


Th? stagnation pressure at the turbine discharge is: 

P56 - PS56 + (W56, PS56, T56) (16) 

The interturbine pressure and temperature are: 


P42 - P56 X fj^p(P4/P5b) 

(17) 

T42 “ T4 X f (P4/P56) 

Ni 

(18) 
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The low pressure turbine exit temperature Is: | 

T56 - TA2 X (P56/P42) (19) 

The total gas flow leaving the mixer and entering the nozzle Is the fan 


flow plus the fuel flow. 

W8 - W2 + WF (20) 

The gas temperature leaving the mixer Is: 

T8 - (W56 X f (T56) + W15 x T15J/W8 (21) 

n 

The nozzle pressure Is: 

P8 - f„ (W8, T8, PO) (22) 

N 

The mixer Inlet static pressure is: 

PS56 - X P8 - X W8^ X T8/P8 (23) 

The fan power Is 

PWl - T2 X W2 X f (P15/P2) (24) 

r* 

and the compressor power is: 

PW25 " C X W21 X (T3 - T21) (25) 

The high and low pressure turbine powers are: 

PW4 - fjjp (T4, P4/P42, W3, WF) (26) 

PW48 - fjj, (T42, P42/P56, W3, WF) (27) 


The above equations describe the static relationships within the en- 
gine. The dynamic states of the engine within the frequency range of the con- 
trols are primarily the rotor accelerations, metering valve and compressor 
variable stator actuators, and the thermal Inertia of the temperature sensors. 
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The rotor accelerations are proportional to the unbalanced power and 
inversly proportional to the polar moment of inertia and speed. 


NL - C_ X (PW48 - PW1)/(J- X NL) (28) 

R L 

NH - C_ X (PW4 - PW25)/(J„ X nH) (29) 

K H 

The metering valve and compressor variable stator actuator velocities are 
proportional to the respective electrical currents from the control logic. 

^ ■ S) " ^ C30) 

ZB - Cj « 16 (31) 


The rate of change of the temperature sensors is proportional to the 
weight flow of the gas past the sensor and the difference between the gas 
temperature and the sensor temperature. 


T3S - X w3 X (T3 - T3S) 


T56S - X W56 X (t56 - T56S) 

TjO 


(32) 

(33) 


The overall accuracy of the model cannot be simply stated because the 
accuracy as operated depends on which states are input. However, the 
accuracy of Individual functions can be given. For example, on the operating 
line the accuracy of the fan and high pressure compressor airflow are with- 
in two percent. Many of the fits are within one percent. 


4.3.2 Feedback Gain and Decision Logic 

One can recognize that the nonlinear model and the feedback-gain matrix 
are actually the well-known extended Kalman filter (Reference 3). Thus the 
feedback gain can be chosen to minimize the mean square error between the 
actual and estimated outputs. The engine can be represented in compact form 
as: 


n + 1 


X + f(x , u ) 
n n n 


g(x , u ) + e' 

® n n n 


(34) 

(35) 
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where, at t • nr, x " state of the engine 


■ input to the engine 

■ neasur.id output 

T 

The neasurraent noise is c* with E(c') • 0 and E(c'c' ) > RI 

n n n n 2 

The "best" estiaate of the neasured output, can be represented by: 


n + 1 


n n n n n n 


- V 


■ “n> - 


(36) 

(37) 


where, at t ■ nr, f^ and represent the nonlinearities in the oodel. 

X “ estinated state of the engine given the past neasurenents 

n 

y^ - expected measured output of the engine 

u ” error in the model which is assumed Gaussian [0,R,] 

n 

c" ■ error in the measurement model which is assumed Gaussian 
" [O.R^'l 

To obtain the best feedback gain, we look at the error, e^^, between 
the states of the engine and the model. 


n*ri n nn isnn n nnn 


n n n n in n n n 


For small errors: 


n + 1 


^e + v - K (y - y) 
nnn •'n ■'n 


y-y»Ce + c 
^n ^n n n 


(38) 

(39) 


(40) 

(41) 


fdtere Gaussian [0, R2I, R2 ■ RA R®i*' K which 

minimizes this error is well known (Reference 3) and is given by the equa- 
tions: 
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( 42 ) 


K - [R, + c p c 

n n 2 n 

^n+l - " K C] P - K C]^ + R, + K R, kI (43) 

n+i n n n 1 n 2 n 


where 


P^ - Efx„ - X) (x„ - i )’') (44) 

n n n n n 

Obviously, the gain matrix given by equation 42 depends on the operating 
point of the engine. 

Under the assumption of Gaussian noise and modeling errors, the error 
between the actual and estimated outputs will also be Gaussian with zero 
mean and covariance. 


R • R, + C P (45) 

y i n 

The expected variance of the error for the i-th Individual sensed variable 

is 


E((y 


i,n 



) (y 


l,n 



T 

) ] 



(46) 


where Ryu Is the 1-th diagonal element of Ry. Thus, the optimum detector 
for failure (Reference 4) is: 




»^R .. 
yii 

No Failure 

(47) 


> k, 

yil 

Failure 

(48) 


The threshold constant Is k. 

When a sensor failure Is detected, the output of that sensor is Ignored 
until it Is again within the tolerance. In the computation of the feedback 
gain, this Is equivalent to setting to zero the corresponding row of the 
matrix C and the row and column of R2 in equations 42 and 43. This Implies 
that the gain K should be computed, at each time step, taking into account 
the available sensors and changes in flight conditions. However, the on- 
board calculation of the Rlccati equation requires considerably more memory 
and computational speed than Is conveniently available. An analysis of the 


33 


•Iganvaluea for th« open and cloaad loopa haa ahown chat a conacanc-fatdback 
utrix can ba found which will yield a atabla ayat«a for the range of power 
at aea level and for varioua eenaor failurea and coabinationa of aenaor 
failurea. The hybrid-coiqNiter aioulation alao confiriM the realiaabilicy of 
of chia aiiq>liflcacion. Thia peraita che K aatrix to be coapuCed off- 
mgine and atill aalncain aatiafactory perforaance. 
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5.0 HYBRID SIMULATION 


The hybrid simulation of the QCSEE OTW model was constructed at the 
General Electric AEG Dynamic Analysis Simulacion Center using two Electronic 
Associates, Inc. (EAI) 690 Hybrid Computing Systems. The engine and basic 
control system were implemented on an EAI 690 Hybrid Computing System which 
consisted of two EAI 680 Analog Computers with a total of 224 amplifiers, an 
EAI 693 Hybrid Interface Unit (Model 1) with 64 channels for A/D (analog/ 
digital) and D/A (digital/analog) conversion, and an EAI 640 Digital Computer 
puter with a 16, 000-word core memory. The FICA portion of the model was im- 
plementei' o i a smaller EAI 690 Hybrid Computing System which consisted of one 
EAI 680 A tog Computoi" v;ith a total of 70 amplifiers, an EAI 693 Hybrid 
Interface Unit ^Model 3) with 32 A/D channels and 24 D/A channels, and an 
EAI 640 Digital Computer with a 15,000-word core memory. Peripheral equip- 
ment includes three EAI 8875 eight-channel strip-chart recorders, two EAI 600 
high-speed line printers, two EAI 500 card readers, an EAI 1700 cartridge 
tape unit, and X-Y recorders. 


5.1 TECHNIQUES 

The two computing systems are linked through trunks between the. analog 
consoles. The computing system containing the engine and basic control 
system is programmed to operate with or without the FICA logic of the second 
computing system. This permitted the design of the basic control system to 
proceed Independently of the FICA development. All engine calculations 
except for noniteration dynamics are performed on the digital computer. 

Bas-'c control calculations are split between tlie digital and the analog 
computers with all dynamics on the analog and most function generation on the 
digital. Transient data is output from the analog computer to the strip- 
chart and X-Y recorders. 

The split of the computational load for the engine and basic control 
system is shown in Figure 13. 

The digital program for simulating the OTW engine and basic control 
system requires approximately 15,300 words of core memory. The resulting 
digital sampling Interval (with steady analog inputs) is approximately 60 
milliseconds. As a consequence, the simulation is run at a time base which 
is twenty times slower than real time. Thus, the digital delay is an effec- 
tive 3 milliseconds in simulation time. 

The FICA portion of the model is primarily a digital simulation requiring 
9,637 words of core memory. The only analog calculations are the temperature 
sensor dynamics and the logic which control the sampling rate and printout 
intervals. The high speed line printer and paper-tape punch are used to 
provide transient output data from the FICA simulation. 
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Figure 13 . Split of QCSEE OTW Simulation 
Computation Load. 



5.2 


S IKULATION VERIFICAT ION 


The simulation of the engine was compared to the cycle deck at several 
operating points which were felt to be of primary Importance for control 
studies. These points were at sea level, static, standard-day conditions 
for the following thrust levels: takeoff, 85% of takeoff, 75% of takeoff, 

65% of takeoff, 50% of takeoff, and 30% of takeoff. The comparison data 
were generated by setting simulation fuel flow and exhaust noc.zle area at 
the cycle deck values. The cycle deck data, simulation data, and percent 
error are showti in Table 1 . 

The simulation data in Table I show good agreement with the cycle 
deck at the power settings from 65 to 100% takeoff thrust. AH percent 
errors are less tlian 2% except for fan stall margin at 65% takeoff thrust. 

At this point, the percent error in fan stall margin is -3.4%; however, since 
the value of stall margin Is 14.91%, the actual difference in terms of 
stall margin is less than 0.5%. Below 65% takeoff thrust, agreement deteri- 
orates to a nuixlmum of 4.81% In not thrust. This deterioration can be at- 
tributed to several factors: round-off errors, truncation errors, and map 
Inaccuracies due to linear Interpolation. The deviation of the model at the 
lower power levels was not considered significant because tin' primary goal 
was to produce a tool which could be used to design a control system which 
would meet the required 1.0-sec accei time from approach to 95% of takeoff 
thrust . 


Since no other source of transient data was available. It was not 
possible to verify the model transiently against an independent data source. 
However, several model transients were run at a lime base one-hundred times 
slower than real time with no observable differences from the 20:1 time 
scale. Thus, the effects due to hybrid simulation (iteration dynamics, A/D 
multiplexer slewing, digital time delay, etc.) were assumed to be negligible. 
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Table 1. Steady-State Verification Data for QCSEE OTW Simulation 
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(Ib/sec) (880.1) (878.8) (820.8) (815.4) 

T8, K 399.9 379.9 -0.02 366.9 367.1 0.06 

(<• R) (683.9) (683.8) (660.4) (660.8) 



Table 1. Steady-State Verification Data for QCSEE OTW Simulation (Continued). 
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6.0 SIMULATION RESULTS 


The QCSEE OTW hybrid simulatiop has been used to develop the control 
system design for the experimental engine. The control system manipulates 
two variables: fuel flow and core compressor stator angle. The system 

Includes a full authority digital electronic control; the design Includes 
three key functions. 

• Reset of the core compressor stator to achieve acceleration- 
thrust-response characteristics required for powerod-lift operation, 

• Fail-fixed servovalve which causes the fuel-metering-valve, power 
piston to f all-in-position if there is a loss in electrical power 
to the digital control or if there is a digital control failure 
resulting in a hard-over, electrical output in cither direction. 

• Failure identification and corrective action logic which recognises 
failures in sensor feedbacks to the digital control and provides 
corrective action for maintaining safe control of the engine. 

Simulation results on each of the above control functions are described 
in the following sections. 


6.1 TRANSIENT THRUST RESPONSE 


The thrust-response requirement for tlie QCSEE Ollv propulsion system at 
altitudes up to 1.8 km (6000 ft) is defined in Figure 14. To simplify tlie dis- 
cussion of transient response, the overall requirement defined in Figure 14 
has been interpreted as a response time from 62 to 95 percent of takeoff 
thrust in 1.0-second maximum. This thrust-response requirement is aimed 
primarily at the landing approach condition where rapid thrust recovery 
capability is required. 

The UTW experimental engine has a variable-pitch fan, variable fan 
bypass duct exhaust nozzle area, and variable core engine compressor stators. 
The variable-pitch fan and variable nozzle area were utilized in tlu' UTK 
control design to .achieve tlie above-cited, tlirust-response requirement. 

However, the OTW experimental engine has a fixed-pitch far. and does not have 
a variable exhaust nozzle area-actuation system. (Note: A variable no;-zJe 

area-actuation system is planned for an OTW flight-type engine.) Early 
simulation studies on the OTW experimental engine predicted thrust-response 
times greater than the 1.0-second requirement. As a consequence, attention 
Wiis directed to more effective utilization of the core compressor stati'rs to 
achieve acceptable thrust-response times with the OTW experimental engine. 

It was determined that the required 1.0-second .acceleration time from 62 to 
95 percent net thrust could be achieved with the following control action. 


Reset the core stators 30 degrees further closed than the normal 
stator schedule at the 62 percent not thrust operating condition. 


dl 




• Rapidly open the core stators to the normal stator schedule during 
the acceleration to 100 percent net thrust. 

Figure 15 Illustrates the potential advantage of using the above core 
stator reset operation. The transient with core stator reset shows that 
acceleration time from 62 to 95% net thrust is less than 1.0 second; whereas, 
the acceleration time with no core stator reset is greater than 1.6 seconds. 


6.1.1 Effect of Core Stator Reset 


Fan speed, and thus engine thrust, during an aircraft approach will be 
controlled by manipulating fuel flow. If the pilot desires the capability 
for fast recovery from approach to takeoff thrust, the proposed procedure is 
to activate the digital control core stator reset mode. The control logic 
causes the core stators to slew In a closed direction. The final stator 
position Is determined by the summation of the core stator reset schedule 
and the normal core stator schedule. (The reset schedule Is a function of 
power setting, and the normal schedule Is a function of core engine corrected 
speed.) Closing the core stators reduces the compressor airflow for a given 
core corrected speed. If fuel flow Is held constant, the engine cycle will 
balance at a lower fan speed. The fan speed control senses the deceleration 
tendency and Increases fuel flow to maintain the scheduled fan speed. As a 
consequence, the core engine rotor settles out at a higher speed when steady- 
state operation Is achieved. The basic objective of the core stator reset 
at the approach power setting Is for the core engine rotor to settle at a 
speed close to the takeoff power setting. Thus, when the pilot demands a 
fast recovery from approach to takeoff thrust, the core engine rotor is 
already near the final speed, and only the fan rotor needs to be accelerated 
to takeoff speed condition. Such operation removes the dynamic effect of 
the core rotor inertia on engine acceleration and achieves fast response 
from approach to takeoff thrust. 

Figure 1.6 shows the effect of core stator reset on pertinent engine 
variables at 62% of takeoff thrust for the sea level, static, standard day, 
zero-bleed condition. When the stator schedule is reset 30 degrees closed, 
core speed Increases to about 13,800 rpm, which Is some 130 rpm above the 
speed at the takeoff power setting. As the stator schedule Is reset from 0 
to 30 degrees closed, fuel flow Increases from 1785 to 2010 kg/hr (3940 to 
4440 Ib/hr), l.e., an increase of 225 kg/hr (500 Ib/hr); specific fuel 
consumption (sfc) Increases from 0.313 to 0.354. In effect, a trade off Is 
being made between the specific fuel consumption and the capability to meet 
the 1.0-second acceleration requirement when 30 degrees of stator reset Is 
used. The simulation data Indicates that this transient response capability 
can be achieved at the cost of 3.77 kg (8.33 lb) more fuel for every minute 
at the 62% thrust approach condition. 

Using 30 degrees of stator reset results in an actual core stator 
position of 24 degrees at 90.4% corrected core speed (see Figure 16). At 
90.4% corrected speed, the nominal stator schedule as a function of corrected 
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core speed calls for the stators to operate on the full-open stop; that 
is, -6 degrees. Thus -6 degrees, plus the 30 degrees reset, results in the 
above -¥2U degree position. 


6.1.2 Simulated Go-Around Maneuver 

The QCSEE control systra is designed to meet the required transient 
response for a go-around maneuver while maintaining safe engine operation. 
Simulation results for an acceleration from 62 to lOOZ net thrust are shown 
in Figure 17. The left-hand portion of Figure 17 shows the specific values 
of selected engine and control variables at the 62Z net thrust condition. 
After the step Increase in power setting, the figure shows: 

e The response time from 62 to 9SZ of takeoff thrust is 0.73 seconds, 
which is well within the 1.0-second requirement. 

e The core stators start opening within the first 0.05 seconds and 
slew to the full-open stop (-6 degrees) at a rate of 60 degrees 
per second, producing a rapid increase in compressor airflow. 

The full-open stop is reached in approximately 0.55 seconds. 

e Fuel flow increases due to the step in corrected fan speed demand 
(scheduled as a function of power setting). During the first 
0.45 seconds, fuel flow is limited by the WF/PS3 acceleration 
fuel schedule. 

e During the interval from 0.45 to 0.85 seconds, the rate of fuel 

flow increase is limited by the calculated core turbine temperature 
control. This control has anticipated that turbine inlet teoq>era- 
ture is approaching, but still below, the reference limit for 
temperature. In this time interval, turbine temperature overshoots 
the steady-state temperature at takeoff by some 122 K (220* R). 
Minimum compressor stall margin is 14.9%. 

e After 0.85 seconds, the fan speed has accelerated to the point 
where it again controls fuel flow. Core speed peaks at 14,200 
rpm. Both rotor speeds, the core stators, and fuel flow settle 
out at the final takeoff operating levels within the next 1.75 
seconds. 


6.1.3 Acceleration Study 

Several factors were evaluated to establish the transient control 
design for the OTW engine. The design process started with the development 
of the WF/PS3 acceleration fuel schedule. This schedule was designed to 
use 5% of the available core compressor stall margin during accelerations 
with no core stator reset. Results from Iritial simulation studies on 
transient core stall margin are plotted versus corrected core speed in 
Figure 18. The figure contains typical stall margin transients for a 


46 


ORIGINAL PAGE IS 
UP POOR QUALHY 








control with no stator reset and for a control with stator reset; satisfactory 
performance is indicated in both cases. 

The rate of removing core stator reset (i.e., the opening race of the 
compressor stators) was also investigated in the initial design studies. 

As shown in Figure 19, rates greater than 60 degrees per second do not 
provide a significant reduction in acceleration time. The adequacy of the 
fuel pump to supply stator actuator flow for the 60 degrees per second rate 
was checked out by the component engineer; his assessment was that this 
rate could be achieved with the current pump design. This 60 degrees per 
second opening rate was chosen for the control design and maintained through- 
out the remainder of the simulation studies. 

Initial design studies also considered how core turbine inlet temperature 
overshoot affects acceleration time from 62 to 95% net thrust. The simulation 
results in Figure 20 show the acceleration time as a function of the delta 
between peak-temperature overshoot and the steady-state temperature of 100% 
net thrust for control designs with and without core stator reset. Based 
on the results in Figure 20, the following conclusions were reached at this 
stage in the control design effort: 

• Unrealistic turbine overtemperatures would be required to meet 
the one-second acceleration requirement when the control design 
does not include core stator reset. 

• Temperature overhsoots of 28 to 33 kelvins (50 to 60® R) will 
produce acceleration times in the 0.9-second range for a nominal 
control design (using 30 degrees stator reset at 62% thrust). 
Therefore, the control design should allow for more temperature 
overshoot to provide more margin with respect to the 1.0-second 
requirement. The judgement was that the above- indicated, 0.1- 
second margin was not sufficient to account for control component 
tolerances. 

The above study was limited to the 33-kelvin (60° R) overtemperature range 
because further refinements to engine simulation scale factors were needed. 
These refinements were accomplished prior to the final simulation studies 
discussed below. 

The OTW experimental engine has a titanium fan; however, a lower 
inertia, composite fan is planned for a flight-type engine. The Initial 
study predictions in Figure 21 indicate an 0.2-second reduction in acceleration 
time when the composite fan is used and when the control design Includes 30 
degrees core stator reset at 62% net thrust. 

In the final acceleration design studies, the core turbine inlet 
temperaturr* control loop was refined so that the fuel control would remain 
in the acceleration fuel schedule mode for a longer period of time during 
the transient. Rate feedback limits were "tuned" to provide the desired 
anticipation for the different control loops. The final simulation predictions 
for the experimental engine acceleration times as a function of initial 
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Figure 19. Transient Response; Accel Time Versus Rate of Removing 
Core Stator Reset. 
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Figure 21 


Initial Predictions of Transient Response for Accel 
Time Versus Core Stator Reset and Fan Rotor Inertia 



thrust level are shovm In Figure 22. The predicted time from 62 to 95% net 
thvust is 0.73 seconds for the nominal control design. 


6.1. A Deceleration Study 

Since there is no transient requirement for decelerations to approach 
power, the primary concern of the deceleration study was safe engine operation. 
The WF/PS3 deceleration schedule was checked out on the model to ensure 
that combustor blowout conditions will not occur. 

A transient for a deceleration from 100 to 62% of takeoff power at sea 
level, static, standard day conditions is shown in Figure 23. For this 
transient, the core stator control logic is set in the "reset on" mode. 

The response time from 100 to 62% thrust is less than 1.0 second. Thrust 
and fan speed both settle -lose to the final values within 2.5 seconds 
after the step decrease in power setting. The core stator control limits 
the rate at which reset is added to the normal core stator schedule. This 
rate limit is approximately 4 degrees per second; thus, it takes approximately 
7.5 seconds to reset the stators 30 degrees closed as shown by the core 
compressor stator error recording in Figure 23. This closing rate limit is 
significantly lower than the 60 degree per second opening rate limit used 
during accelerations. The slow closing rate prevents the simulation from 
operating in a region where the fan hub characteristics are undefined. 


6.2 FAIL-FIXED SERVOV.UVE STUDY 


The purpose of this study was to design the control logic when the 
digital electronic control is connected to a fail-fixed servovalve. The 
hybrid-computer simulation of the engine and control system was used to 
determine the logic. Results from this simulation study indicated that the 
digital control logic should perform the following functions. 

• An inverse function which compensates for both the servovalve 
dead zone and the nonlinear gain of the pulse-width-modulated 
(PWM) driver amplifier. 

• A dynamic method which compensates for servovalve null shift. 

The schematic in Figure 24 describes the location of these functions in 
relation to other digital control logic and also to the driver ampli- 
fier which is connected to the torque motor of the fail-fixed servovalve. 

The left-most block in Figure 24 represents all schedules, dynamics, and 
selection logic for the several controllers which are contained in the 
digital control. The output of this block represents the controller selected 
to regulate engine fuel flow. This output is connected to the blocks 
representing the null-shift compensation and the Inverse function. The 
digital output word from the Inverse function is the signal to the PWM 
driver amplifier, whose output is connected to the fail-fixed servovalve 
assembly. Subsequent paragraphs in this section discuss simulation study 
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Figure 23. Throttle Chop from 100 to 62% Thrust using Linear 
Servovalve and Core Stator Reset. 
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results which ultimately led to the above compensation logic design. Both 
unipolar and bipolar F.M driver amplifiers were considered In this study. 

The Initial study of the fall-fixed servovalve used the characteristics 
of the unipolar-pulse-driver amplifier because it yields better resolution 
from percent digital word to servoflow. An inverse function was developed 
to compensate for the dead zone and to provide characteristics equivalent 
to a linear servovalve from percent digital control word to servovalve 
flow. The control will not perform properly without this compensation 
because of the large dead zone and the use of rate feedback to provide the 
control-loop Integration. For example, at steady-state conditions, the 
rate feedback could be zero and the fan-speed error too small to overcome 
the dead zone. 

The effect of perfect compensation for the dead zone on a throttle 
chop to approach power and a burst back to takeoff power is shown in Figures 
25 and 26. The effect of having the inverse overcompensated for the dead 
zone is shown in Figure 27 where the compensation is 10% greater than the 
dead zone. This produces a step In the digital-control-word/servoflow 
relationship resulting in a very high gain at null and sustained fuel flow 
oscillations at takeoff. Thus, any dead-zone compensation must be sized 
for the minimum expected dead zone. 

When the maximum null shift (±5% of digital output word) for a unipolar 
pulse was used with the dead-zone compensation, the control would not per- 
form approach-power transients. Depending on the polarity of the null 
shift, the control could either decelerate to 62% takeoff thrust or acceler- 
ate back to takeoff thrust, but not both. This was due to the effect of 
the dead-zone compensation which created a dead zone on one side of null 
and a vertical step on the other side of null. This is illustrated in 
Figure 28 for a maximum, positive, null shift. Because of the dead zone, 
the vertical step does not produce an oscillatory system. However, the 
dead zone prevents transients requiring the same polarity of servoflow as 
the dead zone. 

To improve the performance with null shift, the unipolar-pulse-driver 
amplifier was replaced with the bipolar-pulse-driver amplifier which has a 
maximum null shift of ±1% of the digital output word. The effect of no 
dead-zone compensation on a chop to approach power and a burst back to 
takeoff power is shown in Figures 29 and 30. Note the large inaccuracy in 
thrust at approach, 41% instead of 62% of takeoff thrust. This is due to 
the large dead zone and the use of rate feedback to provide the integration 
of the control. Since the system has reached steady state, the rate feedback 
is zero and the error in fan speed is not large enough to overcome the dead 
zone. A similar result is seen when an attempt is made to accelerate back 
to takeoff power with the engine stopping at 62% of takeoff thrust. 

To provide better accuracy with null shift, a compensation method was 
developed which should not appreciably affect transient operation. The 
method uses the difference of expected and sensed metering-valve rates to 
drive a low-gain integrator; the output is added with the digital control 


57 





ORIGINAL PAGE I« 
OF POOR QUALITY 



Stt LtMl, Stittc, 
Sttndird Diy Conditioni 



Ont SKond 


too 

Powtr 

Sitting, 

* 


0 

1.0 

Inlit Itiroit 

Mach 

Nuiniwr 

0 

SO 

Fin Tip 
Stall Margin, 

% 

0 

40 

Con 

Compnssor 
Stator Error, 
digrMt 

•10 

100 

Con 

Cornctid 

Sptad, 

» 

$0 

n 

Con 

Compnstor Inlit 
CorncM Airflow, 
Ib/SK 

30 

IS 

Con 

Coinpretscr 

Pnsjuri 

Ratio 

S 

1.0 

Con 

Compnssor 

Efficlincy 

0.S 



Figure 26. Throttle Buret from 62 to 100% Thruet using Fall-Fixed 
Servovalve with Perfect Dead-Zone Compensation. 
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Figure 27, Takeoff Thruat using Fail-Fixed Servovalve with 
Dead-Zone Compenaation 10% High. 
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Fail-Fixed Servovalve Characteriatica. 
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Figure 29. Throttle Chop to 62% Power Setting using Bipolar-Pulse, 
Fail-Fixed Servovalve with No Dead-Zone Coapensatioit. 
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Figure 30. Throttle Burst to Takeoff Power Setting using Bipolar-Pulse, 
Fail-Fixed Servovalve with No Dead-Zone Compensation. 
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word to counteract the null shift. If the expected rate gain (KER) and 
expected rate lag (xER) exactly match the hardware, there will be no tran- 
sient effect once the null-shift compensation has settled out. The low 
Integrator gain (KDC) reduces the transient effect of the compensation If 
transients are Initiated before It has settled or If the hardware is not 
perfectly matched. If the dynamics used to generate the expected rate are 
well tuned to the hardware, the compensation will respond to a null shift 
as a lag with time constant approximately equal to KFF/ (KER x KOC). See 
block diagrams In Appendix A which contain the above gain- and time-constant 
symbols . 

The specification for the fall-fixed servovalve portion of the control 
was developed from the above results. This specification Includes an 
adjustment for dead-zone compensation; the adjustment (see Appendix A) will 
be used to "tune out" any dead zone during the control system bench tests 
prior to control delivery to the engine. Since the servovalve characteristics 
are syiianetrlc and linear with respect to effective current In the operating 
range, the Inverse function, referred to as f(DIFFS) in Appendix A, will be 
used to linearize effective current with respect to the digital control 
word. This should make the servovalve characteristic linear with respect 
to the digital control word In the operating range. 

The block diagram In Appendix A shows that the digital control logic 
has been designed to connect either to a fall-fixed servovalve or to a 
linear servovalve. The capability to test the digital electronic control 
with each of these two valve types was a control system requirement. Thus, 
when the servovalve type is changed from fail fixed to linear or vice 
versa, a rebuild of the digital control memory Is not required, and testing 
can proceed with little or no delay. 


6.3 FICA SIMULATION RESULTS 


To evaluate the stability and response of failure detectijn and correc- 
tion, it has been simulated on one hybrid computer; an accurate, nonlinear, 
cycle-balance simulation of the QCSEE and nonlinear control system was on a 
second hybrid computer. The failure detection and correction was connected 
between the sensor simulation and the controls simulation. In this manner 
the performance of the system, as it would operate the actual engine, could 
be predicted. 

The development work on the hybrid simulation showed that the model 
could be simplified to eliminate the states representing heat soak. Also, 
the work showed that sensors not otherwise needed could be eliminated; that 
Is, no extra sensors are needed for the system. 

The typical results obtained are shown In Figures 31, 32, and 33 In 
which the power demand to the controls is stepped between 70 percent of 
takeoff and takeoff power. The response shown In Figure 31 Is for the 
system with all sensors functioning and with the failure detection and 
correction functioning. The simulation showed the Kalman filter tracked 
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Figure 32. Computer Trace with PS3 Sensor Pulled (Concluded). 
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the accurate simulation so that at no point in the transient did any sensor 
errors exceed the tolerances. This latter is indicated by the event trace 
in the margin next to the core turbine inlet temperature. When this trace 
steps outward, one or more sensors are indicated as failed. 

The traces in Figures 32, 33, and 34 show PS3, NL, and NH pulled 
(failed) separately. When PS3 is pulled at takeoff power there is no shift 
in the operating point, but there is a change in the transient response. 

The explanation for this is: PS3 is the primary variable for transient 

fuel flow scheduling, but it is not used for steady-state control. With NL 
pulled there is a slight shift in the operating point, and the settling 
transient is changed because this is a controlled variable at takeoff power 
and 70 percent of takeoff power. With NH pulled the speed input to the 
transient fuel scheduling and the compressor variable stator control are 
estimated. As a result there is a small shift in the core stator position. 

The on-en",lne digital control computer has a 12-bit word with double- 
precision capability. The hybrid computer has a 16-blt word with double- 
precision capability, but software is used to mask the least significant 
bits to simulate the digital computer controller 12-bit word. Comparisons 
have been made of the accuracy (as indicated by the shift set point with a 
sensor pulled) and steady-state hunting with a sensor pulled. The results 
indicate that there is a substantial improvement in the hunting with a 16- 
blt computation. This hunting is shown most strongly in the core turbine 
inlet temperature channel as is shown in Figures 33 and 35. The former is 
with a simulated 12-blt control computer and the latter is with a 16-blt 
control computer, and both are with the XNL sensor pulled. These two 
traces Indicate chat there is no significant difference in accuracy or 
stability to disturbances. 

It can be seen from the figures that good engine control is maintained 
while the transient is altered for failed sensors. A fixed-gain feedback 
matrix has been used which shows that it is not necessary to supply a 
separate feedback matrix for each sensor failure and each power range for 
the engine. Although not shown in the computer traces Included herein, 
multiple sensor failures have been demonstrated in which chop and burst 
transients have been controlled adequately. This further demonstrates the 
advantage of the extended Kalman-filter technique in which a fixed-gain- 
feedback matrix is satisfactory for a wide range of conditions. 
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7.0 SUMMARY OF RESULTS 


A hybrid simulation of the OTW experimental engine has been constructed 
and used to develop the control system dyn<unic design. The engine simulation 
was based on the digital cycle deck used to generate the QCSEE preliminary 
technical requirements. Simulation results for throttle bursts from 62 to 
lOOZ takeoff thrust predict that the experimental engine will o»et the 
dynamic requirement of 62 to 95X takeoff thrust in one second. Results for 
transient stall margins, temperature, and inlet Mach number Indicate that 
safe engine operation will be maintained during this engine transient. 
Transient results also predict safe deceleration capability during throttle 
chops. Throttle chops with core stator reset are much slower than without 
reset to ensure safe operation of the fan. 

Simulation results Indicate that operation of the digital control with 
the fall-fixed servovalve can be designed to meet the accuracy and fast- 
response objectives for the OTW experimental engine. 

Simulation of the Failure Identification and Corrective Action (FICA) 
strategy has demonstrated that the strategy built around an extended Kalman 
filter should meet the objectives. Based on the simulation results, the FICA 
will track sea level bursts and chops without giving a false indication of 
sensor failure. Also, based on the simulation results, the FICA will main- 
tain good engine control for sea level bursts and chops with the single 
failure of an engine sensor. 
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DETAILED BLOCK DIAGR AMS AND SCHEDULE/DYNAMIC SPECIFTCATTnws 
for otw digital fuel control 
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digital fuel control which Is to be used on the first build of nw 
experimental engine (Table. II through VII*: 
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Figure 36. QCSEE OTW Digital Fuel Control Block Diagram (Concluded). 















Figure 37. QCSEE OTW Digital Fuel Control Block Diagram for Lag-Hate Feedback with Rate 
Limit in Ml Control and Integrator with Limit for Servovalve Null-Shift 
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Table II. QCSEE OTW Digital Fuel Control Digital Gains 


KDNIT X kNI X kWT X KWTMA 
KNIWR X kWT X KWTMA 

KT41 X kWT X KWTMA 

KTWR X KWT X KWTMA 

KDN2 X KN2M x KWT x KWTMA 

KDN2 X KN2I x KWT x KWTMA 

KN2WR X KWT X KWTMA 

KACC X KWT X KWTMA 

KDEC X KWT X KWTMA 

KDC X kWNB X KWTMA 

KDC (Average) 


Gain 

Per 

0.1225 mA 

rpm 

4.173 mA 

m/sec 

164.3 mA 

In. /sec 

0.2919 uA 

K 

0.5254 uA 


0.8344 mA 

m/sec 

32.85 mA 

In. /sec 

0.2967 uA 

rpm 

(at DN2MXR = 13,200 

rpm) 

0.1618 yA 

rpm 

(at DN2I * 7200 rpm) 

4.173 mA 

m/sec 

164,3 mA 

In. /sec 

10.3 mA 

m 

405.7 mA 

In. 

10.3 mA 

m 

405.7 mA 

In. 

9.86 mA 

cm/sec 

3.881 mA 

In. /sec 

2.66 mA 

cm/ sec 

1.049 mA 

In. /sec 


Table III. QCSEE OTW Digital Fuel Control Digital Gain Distribution 


Note: 


KWT KWTMA 
KDNIT X knI 
KNIWR 
KT41 
KTWR 

KDN2 X KN2M 
KDN2 X KN2I 
KN2WR 

KACC - KDEC 


7.566 

1.619 X 10“^ 

21.709 

6.9437 X io~^ 

4.3417 

3.9215 X 10“^ 
at DN2MXR - 13200 rpm 

2.1384 X 10"^ 
at DN2I - 7200 rpm 

21.709 

53.62 


The above digital-gain distribution is required for the 
lag-rate-feedback limits specified in table of digital 
constants. 








m 




m 


Table IV. QCSEE OTW Digital Fuel Control Digital Scaling. 


Symbol 

Description 

Range 

DDF 

Sensed [PTO - PSll] 

0 to 82.7 kN/m^ 
(0 to 12 psi) 

DPTO 

Sensed PTO 

13.8 to 131.0 kN/m' 
(2 to 19 psia) 

DDP 

DPTO 

Calculated [ (PTO-PSll) /PTO] 

0 to 0.6858 

DT12 

Sensed Fan Inlet Tempera- 
ture 

233 to 344 K 
(420 to 620° R) 

DT3 

Sensed Compressor Discharge 
Temperature 

219 to 861 K 
(395 to 1550° R) 

DT41C 

Calculated HP Turbine 
Inlet Temperature 

256 to 1922 K 
(460 to 3460° R) 

DNIT 

Sensed LP Turbine Speed 

0 to 9595 rpm 

DN2 

Sensed Core Engine Speed 

0 to 15492 rpm 

DPS 3 

Sensed Compressor Discharge 
Pressure 

0 to 2.07 MN/m^ 
(0 to 300 psia) 

DXMV 

Sensed Metering Valve Posi- 
tion 

0 to 2,065 cm 
(0 to 0.813 in.) 


Table V. QCSEE OTW Digital Fuel Control Time Constants. 
Frequency Range From 0.1 to 10.0 Hz 


Symbol 

Description 

Value 

tER 

Inlet Lag In Servovalve Null-Shift 
Compenstation Circuit 

0.075 sec 

tNIWR 

Rate-Feedback Lag in Fan Speed 
Control 

0.50 sec 

TN2WR 

Rate-Feedback Lag in Max Core Speed 
and Idle Core Speed Controls 

0.50 sec 

tTWR 

Rate-Feedback Lag in T41C Control 

0.10 sec 

tMV 

WF RVDT Feedback Demodulator and 
A/D Converter Lag 

£0.01 sec 

tDP 

(PTO-PSll) Sensor, Demodulator, 
and A/D Converter Lag 

£0.02 sec 

tPTO 

PTO Sensor, Demodulator, and A/D 
Converter Lag 

<0.50 sec 

tT12 

Fan Inlet Temperature Sensor, 
Demodulator, and A/D Converter 

£6.00 sec at 48 
kg/sec per m^^C: 
Ib/sec per ft^) 
Airflow Density 

tNIT 

LP Turbine Speed Sensor, Demodu- 
lator, and A/D Converter Lag 

£0.01 sec 

tN2 

N2 Sensor, Demodulator, and A/D 
Converter Lag 

£0.01 sec 

tT3 

T3 Thermocouple, Electrical Sensor 
Demodulator, and A/D Converter Lag 

* 

tPS3 

PS3 Sensor, Demodulator, and A/D 
Converter Lag 

£0.02 sec 

tWTKA 

D/A Converter and Torque Motor 
Driver Amplifier Lag 

£0.01 sec 

tFFSMA 

Fail-Fixed Servovalue Driver 
Amplifier Lag 

£0.01 sec 
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Table V. 


QCSEE OTW Digital Fuel Control Time Constants (Concluded). 



Frequency Range From 0.1 to 

10.0 Hz 

Symbol 

Description 

Value 

tA 

lead Time Constant In Accel Fuel- 
Flow Servoloop 

0.05 sec 

tD 

Lead Time Constant In Decel Fuel- 
Flow Servoloop 

0.05 sec 

tTW 

Lag Time Constant In Calculated WF 
Input to T41C Calculation 

0.05 sec 


* tT 3 = T (Thermocouple) + T(Sensor, Demodulator, and A/D Converter Lag) 
T (Thermocouple) = 0.952(W3/65.29) 

T(Sensor, Demodulator, A/D) s 0.02 sec 
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Table VI. QCSEE OTW Digital Fuel Control Digital Constants. 


Symbol Description Value 

DLN2WR WF Metering Valve Rate-Feedback Limit 0.3175 

in Max and Idle N2 Speed Controls (0.125 

DLTWR WF Metering Valve Rate-Feedback Limit 0.366 

in T41C Control (0.144 


cm/ sec 
in. /sec) 

cm/ sec 
in. /sec) 


Table VII. QCSEE OTW Digital Fuel Control Schedules and Functions. 


1 . 


2 . 


3. 


Corrected Fan Speed Schedule 

NIK - 12125 - 9.3254 x (T12) 

Power Setting (PLA) Schedule 

MPCT - 0.347138 + 0.834489 E-2 x (PU) - 0.181628 E-4 x (PLA)' 


where 

0 < (PLA) < 100% 
Calculated T41 


T41C - X (t3) + K^2 


Where 


(— V' 

\ PS3C j 


245 


■z.mn 


0.226184 E-1 X 


0. 114039 E-3 X 


/ wF 

\ PS3C j 

/wF_y 

I PS3C/ 


( 


-Y 


WF 

ps3cy 


4. Calculated Engine Fuel Flow as a Function of Calculated Metering- 
Valve Position (XMV) 

For 0 (XMV) <0.2809 

WF - 162.2 + 1653.33 x (XMV) + 18277.8 x (XMV)^ 


For 0.2809 < (XMV) £ 0.813 

WF - -348.68 + 4151.59 x (XMV) + 15864.5 x (XMV)^ 

Note: These Equations are based on the average of test data from WGC 

Hydromechanical controls S/N 12359 and S/N 13082; the following 
equation defines the relation between metering-valve position 
(XMV) and RVDT feedback volts/volt 


H9 


Table VII. QCSEE OTW Digital Fuel Control Schadulea and Functions 
(Continued) . 


Calculated /ei2 

For 244 < T12 < 322 K (440 < T12 < 580’ R) 

RTH12 - 0.383323 + 0.731799 x (TH12) - 0.115215 (TH12)^ 

Max Corrected Core Speed Schedule 
For 0 (PLA) < 46.57Z 

N2K - 7200 + 153.53 x (PLA) 

For 46.57 < (PLA) < 100% 

N2K - 14350. 

m*y g ^ 

Calculated ' T 12 " " ** * function of Percent Corrected Fan Speed 
For 0 < PCNIR < 39.862 

2 

DTQT12 - 0.358845 E-2 + 0.195553 E-3 x (PCNIR) + 0.547267 E-5 x (PCNIR) 
For 39.862 < PCNIR <62.8 

2 

DTQT12 • 0.895999 E-2 - 0.518285 E-3 x (PCNIR) + 0.200000 E-4 x (PCNIR) 
For 62.8 < PCNIR < 94.13 

DTQT12 - -0.249999 E-2 + 0.559521 E-4 x (PCNIR) 

2 

+ 0.137619 E-4 X (PCNIR) 

For 94.13 i PCNIR 

DTQT12 - 1.3351 - 0.273429 E-1 x (pcNlR) + 

0.15387 E-3 X (PCNIR)^ 

Calculated /625 

For 244 <. T25 < 379 K (440 < T25 < 683* R) 

RTH25 - 0.383323 + 0.731799 x (TH25) - 0.115215 x (TH25)^ 


Table VII. QCSEE OTW Digital Fuel Control Schedules and Functions 
(Continued) . 


9. Ar.cel Fuel Schedule 
If XN < 2600 

PHI - 15.96 
If XN i 2600 

PHIBDY - XN fA2 + XN X (A^ + A^ x XN)] 

PHIT41 - 69.03 - 0.065 * T25 - 0.0017 (XN - 12700) 
PHI ■ AMINl (PHIBDY, PHIT41) 


XN 


^2 

^3 

^4 

2,600<XN<6,800 

O.42118701E2 

-0.25707655E-1 

0. 74804485E-5 

-0.56754241E-9 

6,800£XN<10,20 

-0.19708559E3 

0.86346402E-1 

-0.10384008E-4 

0.39665754E-9 

10,200£XN 

0.36912564E3 

-0.80335893E-1 

0.59159283E-5 

-0.13263709E-9 


XN • Core rpn//T25/288.16 (rpm/»/T25/518.69) 

T25 * Core Inlet Temperature, K (* R) 

PHI - WF/[^PS3C X (125/518.69)°*^®®] 

10. Calculated (02j)°*®®® 

For 266 < T25 < 379 K 
(478 £ T25 < 683* R) 

TH688 - 0.211270 + 0.888286 x (TH25) - 0.996309 E-1 x (TH25)^ 

11, 12, 13. Conversion of scheduled accel WF, decel WF, and sensed WF to 

equivalent WF metering valve position. These three functions Must Use 
the same polynomial -equation fit, which is: 

XMV - 0.31297E-1 + 0.143045E-3 x (WF) - 0.139539E-7 x (WF)^ 

+ Q.613485E-12 x (WF)^ 
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Table VII. QCSEE OTW Digital Fuel Control Schedulaa and Functions 
(Concluded) . 


This equation is based on average of test data from WGC hydro- 
mechanical controls S/N 12359 and S/N 13082 and on the WGC 
equation defining the relation between metering valve position 
(XMV) RVDT feedback volts/volt, i.e: 

(XMV) - 0.813 + 


1 


92 


APPENDIX B 


DETAILED BLOCK DIAGRAM AND SCHEDULE /DYNAMIC SPECIFICATIONS FOR OTW DIGITAL 

CORE STATOR CONTROL 


Includt-d in this appendix are the block diagram (Figure 39} and the 
specifications for gains, time constants, schedules, and limits (Tables VIII 
through XII) for the digital core stator control which Is to be used on the 
first build of the OTW experimental engine. 
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Schedule 











Table VIII. OTW Core Stator Digital Control Digital Scaling. 


Symbol 

Description 

Range 

DZB 

Core Stator Position Feedback 

-0.137 to +4.895 cm 
(-0.054 to +1.927 in.) 


Table IX. OTW Core Stator Digital Control Time Constants. 
Frequency Range From 0.1 to 10.0 Hz 


Symbol 

Description 

Value 

tZTMA 

D/A Converter and Torque Motor 

£ 0.01 sec 


Driver Amplifier Lag 


tZB 

LVDT Demodulator and A/D 

< 0.01 sec 


Converter Lag 



Table X. OTW Core Stator Digital Control Digital Constants. 


Symbol 

Description 

Value 

DZBROR 

Core Stator Reset Opening- 

-2.79 cm/sec 


Rate Limit 

(-1.10 in. /sec) 

DZBRCR 

Core Stator Reset Closing- 

+0.19 cm/sec 


Rate Limit 

(+0.075 in. /sec) 
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Table XI. QCSEE OTW Core Stator Schedule 


e FB - fi 4 (XN) 

5 Second-Degree Equations 


10 DIMENSION A(3.5) 

20 DATA A 

30&/Z.2795249I - .0131*51068. . 71650476 E- 4, 
408,.87748352» .0326 17459, -.30322 192E-3, 
50&8. 3086021 I- . 14282437, .72784362E- 3, 
608.-5.7799969, . 16498718, -. 952 1339 lE-3 , 
708.16.450925,- .22620007, .76669457E- 3 
308/ 

90 900 FORMAT (V) 

100 DO 100 Jsl,66 

110 READ(1S,900)L,XN,FBACT 

120 1=5 

130 IF(117,-XN) 10,10,1 
140 L 1=4 

150 IF(94.-XN) 10, 10,2 
160 £ 1=3 

170 IF (79. -XN) 10, 10,3 
180 3 1=2 

190 IF(61 .-XN) 10, 10,4 

200 4 1=1 

210 10. FB = A(1 , I )+XN* <A(2, I ) +A(3, I)*XN) 
220 FB=AMAXI (0. ,AM1N1 (FB, 1 .8688) ) 

230 FBERR=FBACT-FB 

240 100 WRITE(14,900) XN,FB,FBERR 

250 ST0P5END 


Coefficients 

4 A(1,1).A(2,1),A(3,1) 
4 A(1,2),A(2,2).A(3,2) 
4 A(1,3),A(2,3),A(3,3) 
4 A(1,4).A(2,4).A(3,4) 
4 A(1,5).A(2,5),A(3,5) 


, Coefficient Selection 
Log 


4 Feedback Equation 
4 Min & Max Schedule Flats 


NOMENCLATURE 


FB “ Feedback Position in inches 


XN - 


Core rpm 



418.69 


X 10-2 



Note 
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Table XII. QCSEE OTW Core Stator Reset Schedule. 


DELZB - 

(% PLA) 


For 0 < PLA s 80% 

DFXZB -0.77 


For 80 < PLA i 90% 

DELZB = -0.77E-1 

X PLA + 6.93 

For 90 < PLA i 100% 

DELZB =0.0 



APPENDIX C 


A8 

AES 

ALT 

BETA 

C 

CPI 3 
CP3 

DLBETA 

DTAMB 

FAR41 

FAR49 

FAR5 

FAR6 

FN 

GAM6 

HI 2 

HI 3 

H3 

H36 

HAl 

H49 

H5 

H6 

HCL42 

HCU8 

HCL49 

HCL5 

HCL55 

IB, IS 

IWF, IWT 

K 

KC 

KF 


N2, NH 
NL 

NHPLA 

PO 

Pll 

P12 

P13 

P15 


NCaiENCLATURE 


Exhaust nozzle throat actual area 
Exhaust nozzle throat effective area 
Altitude 

HP compressor stator angle 

Various subscripted constants 

Fan tip discharge specific heat 

HP compressor discharge specific heat 

HP compressor stator error 

Adder on standard atmosphere temperature 

HP turbine rotor Inlet fuel/alr ratio 

LP turbine rotor Inlet fuel/alr ratio 

LP turbine discharge fuel/alr ratio 

Mixer discharge fuel/alr ratio 

Net thrust 

Mixer discharge specific heat ratio 

Fan tip inlet enthalpy 

Fan tip discharge enthalpy 

HP compressor discharge enthalpy 

Combustor inlet enthalpy 

HP turbine rotor Inlet enthalpy 

LP turbine rotor Inlet enthalpy 

LP turbine discharge enthalpy 

Mixer discharge enthalpy 

Enthalpy of WACL42 

Enthalpy of WACL48 

Enthalpy of WACL49 

Enthalpy of WACL5 

Enthalpy of WCL55 

Torque motor current from digital stator 
control 

Torque motor current from digital fuel 

control 

Gain 

Gain proportional to reciprocal of core 
rotor moment of Inertia 
Gain proportional to reciprocal of total 
rotor moment of Inertia for fan, gearbox, 
and fan turbine 

HP compressor rotor physical speed 

Fan physical speed 

Core speed demand from H/M PLA cam 

Free-streara pressure 

Inlet throat total pressure 

Fan tip Inlet total prei’sure 

Fan tip discharge total pressure 

Bypass duct inlet total pressure 


m2 (in. 2) 
m2 (in. 2) 

m (ft) 
degrees 

J/kg-K (Btu/lb-“ R) 
J/kg-K (Btu/lb-“ R) 
degrees 
K (“ R) 


N (Ibf) 

J/kg (Btu/lb) 
J/kg (Btu/lb) 
J/kg (Btu/lb) 
J/kg (Btu/lb) 
J/kg (Btu/lb) 
J/kg (Btu/lb) 
J/kg (Btu/lb) 
J/kg (Btu/lb) 
J/kg (Btu/lb) 
J/kg (Btu/lb) 
J/kg (Btu/lb) 
J/kg (Btu/lb) 
J/kg (Btu/lb) 
mA 


mA 


rpm/sec/kg-m 
(rpm/sec/lb-f t) 
rpm/sec/kg-m 
(rpm/sec/lb-ft) 

rpm 

rpm 

rpm 

N/m2 (psla) 

N/m2 (psla) 

N/m2 (psla) 

N/m2 (psla) 

N/m2 (psla) 


OF POOH Ol’AniT^ 


98 


PI 6 

Bypass duct discharge total pressure 

N/m2 (psia) 


P2 

Fan hub inlet total pressure 

N/m2 (psia) 


P21 

Fan hub discharge total pressure 

N/m2 (psia) 


P23Q2 

P25 

Fan hub pressure ratio 
HP compressor inlet total pressure 

N/m2 (psia) 


P3 

HP compressor discharge total pressure 

N/m2 (psia) 


P4 

Combustor discharge total pressure 

N/m2 (psia) 


PA 2 

HP *-urbine discharge total pressure 

N/m2 (psia) 


PAS 

LP turbine vane inlet total pressure 

N/m2 (psia) 


P5 

LP turbine discharge total pressure 

N/m^ (psia) 


P56 

Core duct discharge total pressure 

N/m2 (psia) 


P6 

Mixer discharge total pressure 

N/m^ (psia) 


P8 

Exhaust nozzle throat total pressure 

N/m2 (psia) 


PCNA 

PCNLR 

PHI12 

Percent HP compressor physical speed 
Percent fan corrected speed 
Fan tip inlet entropy function 

J/kg-K (Btu/lb-° 

R) 

PHIAl 

HP turbine rotor inlet entropy function 

J/kg-K (Btu/lb-° 

R) 

PHIA9 

LP turbine rotor inlet entropy function 

J/kg-R (Btu/lb-° 

R) 

PLA 

PSll 

Power lever angle (power setting) 
Inlet throat static pressure 

N/ni2 (psia) 


PS16 

Bypass duct discharge static pressure 

N/m^ (psia) 


PS3 

HP compressor discharge static pressure 

N/m2 (psia) 


PS3MC 

Control-sensed PS3 

N/m2 (psia) 


PS3MEC 

Control PS3 after FICA logic 

N/ra2 (psia) 


PS56 

Core duct discharge static pressure 

N/m2 (psia) 


PS8 

Exhaust nozzle throat static pressure 

N/m2 (psia) 


PWl 

PWl 2 + PW2 

W (hp) 


PW12 

Fan tip power 

W (hp) 


PW2 

Fan hub power 

W (hp) 


PW25 

HP compressor power 

W (hp) 


PWA, PWAl 

HP turbine power 

W (hp) 


PWA8, PWA 9 

LP turbine power 

W (hp) 


PWPXH 

HP rotor power loss 

W (hp) 


PWPXL 

LP rotor power loss 

w (hp) 


R12 

Fan tip inlet gas constant 

m/K (ft/“ R) 


S 

Laplace variable 

sec~l 


SLS 
SMI 2 
SM2 
SM25 

Sea level, static condition 
Percent fan tip stall margin % constant flow 
Percent fan hub stall margin constant flow 

Percent HP compressor stal 1 margin % constant 


TO 

f low 

Free-stream temperature 

K C R) 


Til 

Inlet throat total temperature 

K (° R) 


T12 

Fan tip inlet total temperature 

K C R) 


T13 

Fan tip discharge total temperature 

K (“ R) 


T15 

Bypass duct inlet total temperature 

K r R) 


T2 

Fan hub inlet total temperature 

K C R) 


T21 

Fan liub discharge total temperature 

K (“ R) 


T25 

HP compressor inlet total temperature 

K C R) 


T3 

HP compressor dlscljarge total temperture 

K C R) 


T36 

Combustor inlet total temperature 

K (“ R) 
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T35 

Control-sensed T3 

K (• R) 

T4 

Combustor discharge total temperature 

K (• R) 

T41 

HP turbine rotor inlet total temperature 

K <• R) 

TAIC 

Control-calculated T41 

K (• R) 

T42 

HP turbine discharge total temperature 

K (• R) 

T48 

LP turbine vane inlet total temperature 

K (• R) 

T49 

LP turbine rotor inlet total temperature 

K (• R) 

T55 

LP turbine frame discharge total 
temperature 

K (• R) 

T56 

Core duct discharge total temperature 

K (• R) 

T56S 

Control-sensed TS6 

K (• R) 

T6 

Mixer discharge total temperature 

K (" R) 

T8 

Exhaust nozzle throat total temperature 

K (• R) 

W12 

Fan tip inlet airflow 

kg/sec (Ib/sec) 

W15 

Bypass duct inlet airflow 

kg/sec (Ib/sec) 

W2 

Fan hub inlet airflow 

kg/sec (Ib/sec) 

W21 

Fan hub discharge airflow 

kg/sec (Ib/sec) 

W25 

HP compressor inlet airflow 

kg/sec (Ib/sec) 

W2A 

Fan front-face total flow 

kg/sec (Ib/sec) 

W3 

HP compressor discharge flow 

kg/sec (Ib/sec) 

W36 

Combustor inlet gas flow 

kg/sec (Ib/sec) 

W4 

Combustor discharge gas flow 

kg/sec (Ib/sec) 

W41 

HP turbine rotor inlet gas flow 

kg/sec (Ib/sec) 

W48 

LP turbine vane inlet gas flow 

kg/sec (Ib/sec) 

W49 

LP turbine rotor inlet gas flow 

kg/sec (Ib/sec) 

US 

LP turbine discharge gas flow 

kg/sec (Ib/sec) 

use 

Core duct discharge gas flow 

kg/sec (Ib/sec) 

U6 

Mixer discharge gas flow 

kg/sec (Ib/sec) 

U8 

Exhaust nozzle throat gas flow 

kg/sec (Ib/sec) 

UA41 

HP turbine rotor inlet airflow 

kg/sec (Ib/sec) 

UA49 

LP turbine rotor inlet airflow 

kg/sec (Ib/sec) 

UAS 

LP turbine discharge airflow 

kg/sec (Ib/sec) 

UACL42 

HP turbine discharge cooling flow 

kg/sec (Ib/sec) 

UACL48 

LP turbine vane inlet cooling flow 

kg/sec (Ib/sec) 

UACL49 

LP turbine rotor inlet cooling flow 

kg/sec (Ib/sec) 

UACLS 

LP turbine discharge cooling flow 

kg/sec (Ib/sec) 

UCLSS 

LP turbine frame discharge cooling flow 

kg/sec (Ib/sec) 

UF, UFM 

Engine fuel flow 

kg/hr (Ib/hr) 

XM 

Flight Mach munber 


XMll 

Inlet throat Mach number 


XMV, ZwF 

Fuel metering valve power piston stroke 

cm (in.) 

XNH 

HP compressor physical speed 

rpm 

XML 

Fan physical speed 

rpm 

ZBETA, ZB 

HP compressor stator actuator stroke 

cm (in.) 

z-1 

Z transform 


T 

Time constant 

sec 
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